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Abstract. We present SO and SO 2 observations in the region of the low mass 
protostar IRAS2/NGC1333. The East-West outflow originating from this source 
has been mapped in four transitions of both SO and SO 2 . In addition, CS obser¬ 
vations published in the literature have been used. We compute the SO, SO 2 and 
CS column densities and the physical conditions at several positions of the outflow 
using LTE and a non-LTE LVG approximations. The SO 2 /SO and CS/SO abun¬ 
dance ratios are compared with the theoretical predictions of a chemical model 
adapted to the physical conditions in the IRAS2 outflow. 

The SO 2 /SO abundance ratios are constant in the two lobes of the outflow 
whereas CS/SO is up to 6 times lower in the shocked gas of the East lobe than 
in the West one. The comparison with the chemical model allows us to constrain 
the age of the outflow produced by IRAS2 to < 5 x 10® yr. We find low densities 
and temperatures for the outflow of IRAS2 (< 10® cm“® and < 70 K) from SO 
and SO 2 emission probably because the two molecules trace the cooled entrained 
material. The East lobe of the outflow shows denser gas compared to the West 
lobe. We also discuss some constraints on the depleted form of sulphur. 
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1. Introduction 


Accretion and ejection are two apparently antithetic phenomena occurring simultane¬ 
ously during the first stages of star formation. Both are vital in the overall process: 
the former for the central object to grow, the latter to eliminate the excess angular 
momentum, and to allow accretion to proceed. In the very first phases of star forma¬ 
tion, the ejection occurs through collimated, fast and extended flows of material, of¬ 
ten probed by molecular rotational lines, and therefore refered to as molecular outflows 
and/or jets. Indeed, the forming star being still embedded in the parental molecular 
cloud, the outflowing gas strikes the ambient material, often causing molecular shocks 
at the interaction interface. These shocks can be and indeed often are so violent that 
the grain mantles are partially destroyed, so that some heavy elements, usually frozen 
onto the grain surfaces, are released into the gas phase. Depending on the shock velocity, 
magnetic field and pre-shock density, the s hocks may dissociate the existing molecules 
(J-type shocks: iHollenbach fc McKeelll989^ . or, on the contrary, form many molecules 


in th e “gently” shocked gas (C-type shocks 


: lDrainj[l^ : 


Flower fc Pinea.ii des Foreti 


200,11) . or in the post-shocked g as. The SiO molecule represents a representative case in 
this respect i Schilke et al .111 997 ). but several other molecules have been observed to have 


abundances much greater in molecular shocks than in molecular clouds (e.g. 


Bachillei 


199m . The time scale of the ejection process compared to the protostar life is not well 
known. To date the outflows there are two possible methods. One is dynamical and uses 
the ejected matter velocity (based on the line profiles) and the distance from the pro¬ 
tostar. An other way to constrain the age of shocks is to study their chemical evolution 
through chemical clocks. Chemical clocks are abundant species whose abundance depends 
more on time compared with the dynamical evolution of the source. Thus comparing the 
observed abundances of these species with their theoretical evolution should indicate the 
age of the shocked region. 

Sulphur-bearing species are of particular interest, for the following two reasons: 
a) Sulphur is known to be severely depleted in molecular clouds, where its me asured 


abundance is 1000 times less than the cosmic abundance l|Tieftriink et al 


1994 . Note 


that the sulphur depletion likely occurs in the molecular cloud phase because in the 
diffuse medium the total abundance of sulphur in the gas phase is close to the cosmic 
abundance of S. On the other hand, SO, SO9 and H^S are often abundant in molecu¬ 


lar shocks llBachiller fc Perez Gutierrez 


1997 


Codella fc Bachillei 


199£; 


Wakelam et al 


2OO41J) : therefore S-bearing molecules can be potentially good shock tracers; 
b) Sulphur chemistry is relatively fast, with typical timescales of some 10"'’ yr 


Wakelam et al 


2004a|) : therefore, the relative abundanc es of “appropriate” S -bearing 


molecules could in principle be used to date the outflows llBachiller et al 


2001 ). 
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have focused on the use of sulphur-bearing species as “chemical clocks” J 

Charnlev 

1997 

Hatchell et al. 

1995 

Bachiller et al. 

2001, 

van der Tak et al. 

2003 

Wakelam et al. 

2004blall. Several observational studies have targeted the “hot cores” of protostars, 


regions where the dust mantles evap orate because of the elevated dust temperature. 


Unfort unat ely, 


Wakelam et a,1 


2flfl4ari have demonstrated that the use of S-bearing mo¬ 


lecules as chemical clocks is far from easy. The main difficulty is that, indeed, the time 
evolution of the S-bearing species abundances critically depends on the initial form of 
the sulphur on the grain mantles, and on the exact physical conditions of the gas. As a 
result, it is very difficult to estimate the age of sources belonging to different molecular 
clouds since the initial sulphur composition could be very different and mask any evolu¬ 
tionary effect on the S-bearing species abundances. It is, however, still possible that the 
abundance of S-bearing molecules in shocked gas belonging to the same system (outflow 
source and molecular cloud) can be used to date the different shock sites. At least, one 
can expect that this would be easier to interpret than in the case of hot cores of different 
sources. In other words, a system where ejection events occurred at different times and 
caused successive shock sites could be a good benchmark for a time-dependent sulphur 
chemistry study. For this reason we decided to study an outflow system where multiple 
ejection events are suspected to have occurred. The basic question we want to answer is: 
can the (relative) abundance of sulphur-bearing molecules identify shocks that occurred at 
different times? The article is organized as follows. We describe the selected outflow sys¬ 
tem (§2), the observations (§3), the obtained maps (§4) of the most abundant S-bearing 
molecules in the region and the de rived relative abundan ces along the outflow (§5). In 
§6 we describe the chemical model I Wakelam et ahlLnOda ) applied to the outflow system 
studied. Finally, in §7, we discuss the results and compare the observed abundances with 
the predictions of the chemical model. 


2. NGC1333-IRAS2: source background 


The selected source, IRA S2, is situated in the NGC1333 molecular cloud, an active star 


forming region at 220 oc 

jCernif 

1990 

)i "w 

stars ( 

Asoin et al. 

1994 

; 1 

_/ada et al. 

1997) 


IflflT ). with several low- and intermediate-mass proto- 


1997ll . Evidence of the interaction between the pro- 


have been widely reported ( 

Bachiller Kr. CermcharnI 

1990 

Sandell et alll994: 

Bla.ke et al. 

1995; 

Wariu et al. 

1996; 

banger et al. 


Va,r(i-ThoTrir)Son et a,l. 

1996 

Ijefloch et al. 

1998a 


Bachiller et al. 

1995 

Knee & Sandel 

hoof 

1 ). NGC1333 hosts a few very young 


Class 0 sources only detectable in the millimeter to FIR wavelength range ( .lennings et al 
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Among them, IRAS2 (IRAS 03258+3104 source of the IR AS Point Source Ca talogue) 
is a low mass Class 0 protostar located at the edge of a cavity ( Langer et alJlQQO l. IRAS 2 
is indeed composed of three continuum components: A, B and C. The B and C compo¬ 
nents (IRAS2B and IRAS2C) a re respectivelly at 30" South-East and 30 " North-West 


of the A component (IRAS2A) IILoonev et al 


200C; 


Knee fc Sandell 


200(11 . In practice, 


IRAS2A is the most studied source among the three, and it is often refered to as IRAS2 


I Sandell et al 


1 904jl . In the following, we will follow this tradition of calling IRAS2 the 


A component. 


Two bipolar outflows in the North-South and East-West directions emanate from IRAS2, 

suggesting that it is a non-resolved binary system ( 

Sandell et al. 

1994 

). The North- 

South outflow was first detected by 

Liseau et al. 

(1988 

) in the CO 1 

- 0 transition. 


Later on. 

Sandell et al. 

Knee fc SandelJ 

2001 ] 

, 


1 199411 showed that this outflow is not very collimated (see also 
for a detained study of this outflow). These authors also discov¬ 
ered another outflow coming out from IRAS2 and oriented East-West. This outflow is 


the East and West lobes spatially well separated I 

Bachiller et al. 

1998 

). Using interfere- 

metric maps of the CH 3 OH emission. 

Bachiller et al. 

1 1998ll have proposed that the lobes 


of this outflow are composed of several bullets, created in successive mass loss episodes. 
They predicted a difference of ~ 2 x 10^ yr between the farthest and nearest bullets. 

We obtained maps of the IRAS2 East-West outflow system in four transitions of SO 
and SO 2 respectively. The first goal is to estimate the gas temperature and density as 
well as the column density of SO and SO 2 along the outflow. We also re-computed the 
C S column dens i ties a long the outflow, using the line emission maps previously published 
1 199(ir . The ultimate goal is to compare the SO /SO 2 /CS ratios observe d 


banger et a,l 


along the outflow with the predictions of the chemical model by 
and to estimate the age of the shocks along the outflow. 


Wakelam et al 


(2004a|), 


3. Observations 

The IRAS2 region was observed with the IRAM 30m telescope at Pico Veleta (near 
Granada, Spain) during several runs in January and August 1998 and January 1999. We 
performed large-scale maps at a sampling of 12 " in the following molecular transitions: 
SO (32 —> 2 i, 23 —> I 2 , 43 ^ 32 and 65 —> 64 transitions) and SO 2 ( 3 i ^3 ^ 2 o, 2 , 
10 i ,9 ^ lOo.io, 5 i ^5 ^ 4 o ,4 and 62,4 ^ 4 i ^3 transitions). We covered an area of Aa x AJ 
= 228"X 96" in SO and 216"x 84" in SO 2 . The relative coordinates (Aa, AS) used in all 
the maps presented here refer to the source IRAS2 (a(2000.0) = 3’^28™55®41, J(2000.0) 
= 3l'’l4'35.08"). TableHsummarizes the relevant parameters of the observed transitions. 
All observations were performed using the position switching mode, with an off posi- 
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Table 1. Parameters of the observed SO and SO 2 transitions: the frequency {i/), the 
energy of the upper level divided by the Boltzmann constant (E„p/k), the telescope 
beam, the beam efficiency of the telescope (Beff) and the spectral resolution (Ai/). 


Transitions 

u 

(GHz) 

Eup /k 

(K) 

Beam 

// 

Beff 

Al/ 

(kHz) 

SO 32 ^ 2i 

99.299 

9 

24 

0.73 

40 

SO 23 ^ I 2 

109.252 

21 

22 

0.68 

40 

SO 43 —^ 82 

138.178 

16 

18 

0.54 

40 

so 65 ^ 54 

219.949 

35 

11 

0.41 

80 

SO 2 3i,3 ^ 2o,2 

104.029 

8 

24 

0.69 

40 

SO 2 10 l ,9 ^ lOo.lO 

104.239 

55 

24 

0.69 

40 

SO 2 5i,5 ^ 4o,4 

135.696 

16 

19 

0.56 

40 

SO 2 62,4 ^ 4i^3 

241.615 

24 

11 

0.39 

80 


This position is void of molecular emission si nce it corresponds to the center of the molec- 

( 199fir . Pointing was checked every 1-2 hrs on the 


ular cavity described in 


Lang-er et al. 


nearby continuum source 0333-1-321. The pointing accuracy was found to be better than 
S". An autocorrelator split into three parts was connected to three receivers working at 
1, 2 and 3 mm respectively, allowing simultaneous observations of three different transi¬ 
tions. The achieved spectral resolution is < 0.12 km s“^. The weather conditions were 
very good, and typical system temperatures were in the range of 120, 200 and 300 K at 3, 
2 and 1.3 mm respectively. The spectra were calibrated with the standard chopper-wheel 
method and are reported here in units of main-beam brightness temperatures. 


4. Results 


Figured displays the integrated intensity maps of the SO 32 ^ 2i and 43 ^ 32 lines, as 
well as the S O 9 3i .3 ^ 2n 9 and 5 i ^ 4o,4 lines, and, for comparison, the CS 2 ^ 1 and 


2 lines llLanger et al. 


199^. Three peaks of emission stand out on the maps: one 


associated with the East lobe of the outflow (marked E in the map - offset 72", -12") and 
two associated with the West lobe (marked W1 -offsets -60", 12"- and W2 -offset -96", 
24"- respectively). On the contrary, no enhanced emission is observed towards the source 
itself with the exception of the CS 3^2 line. The E and W2 r egions are symmetr ically 

1 1996 1 ^ as 


situated around IRAS2, and have previously been identified by 


Langer et al 


the bow shocks at the end of the East and West lobes of t he outflow. The three positions 


E, W1 and W2 coincide with peaks of CH 3 OH emission llBachiller et al 


1998). 


For a proper comparison of the spectra taken at different frequencies, all spectra were 
smoothed to the same spatial resolution of 24", which corresponds to the largest beam 


^ The CS 5 —> 4 line has also been observed but is seldon detected. Thus, we did not use it 
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Aa [arcsec] 

Fig. 1. Integrated intensity maps of the CS 2 ^ 1 (a) CS 3 ^ 2 (b) SO 82 ^ 2i (c) SO 
43 ^ 82 (d) SO 2 3 i ,3 —> 2 o ,2 (e) and SO 2 5 i _5 ^ 4 o ,4 (f) transitions, (a), (b), (c) and 
(d) : first level is 0.5 K km s“^ with level step of 2.5 K km s“^. (e) and (f ): first level is 
0.2 K km with level step of 0.3 K km s“^. The CS data are taken from lT;a,nger et a,l 


(122^ whereas the SO and SO 2 data are from this work. In each map, the black points 
represent the observed positions. The star shows the position of the protostellar source 
IRAS2. The arrows point to the positions whose spectra are displayed in Fig. |21 


(including the CS observations from 


Langer et al 


199111 . Note that the SO 2 5 i _5 ^ 4 o ,4 


and 52,4 —> 4 i ,3 spectra at the IRAS2 position were not smoothed because of the lack 
of observed positions around the source. However, since the target of the present study 
is the outflow and not the source itself, this is not relevant. The spectra of the SO and 
SO 2 transitions observed towards IRAS2, E, W1 and W2 are shown in Fig. 12 and the 
derived line parameters are summarized in Table 12 

As already noticed, the emission of the SO and SO 2 lines is strongly enhanced in the E, 
W1 and W2 positions, whereas it is very weak toward the protostar itself. In the direction 
of IRAS2, the spectra exhibit profiles similar to the ones observed in the molecular cloud: 
very narrow lines with an intensity decreasing with the energy of the line suggesting that 
the emission is dominated by the cloud rather than the source. Conversely, the spectra 
of the outflow regions are very intense and exhibit high-velocity wings. In the western 
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(—60", 12") and W2 (—96", 24") (see text). All shown spectra are smoothed at the 
largest beam 24". The SO 2 5 i ,5 —> 4 o ,4 and 62,4 ^ 4 i _3 toward IRAS2 are not repre¬ 
sented because they cannot be smoothed (see text). The vertical dashed line on each 
spectrum marks the position of the cloud systemic velocity (vlsr = 7.6 km s“^). 


component, with a broad profile, is blue-shifted by 6 km s“^ with respect to the velocity 
of the ambient cloud. We will refer to this component as the High Velocity Component, 
“HVC”. The second component peaks at the systemic velocity of the cloud, 7.6 km s“^. 
We will refer to this component as the Low Velocity Component, “LVC”. While the 
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Fig. 3. Integrated intensity maps of the LVC and HVC SO da ^ 82 transition separately 
traced on the left and right panels respectively. The first level is 0.5 km s“^ with level 
step of 1 km s“^. 


component because the maps are not large enough. However, there is some indication 
that the LVC component is associated with the molecular cloud, for the LVC component 
is extended and is poorly contrasted (see Fig. El). Also, the vlsr, 7.6 km s is constant 
across the emitting region. The lines of the East lobe (E) are brighter than those of the 
West lobe, and only one component is detected, red-shifted with respect to the systemic 
velocity. Most probably the component centered on the cloud velocity is buried of this 
red-shifted bright component, making it difficult to separate the two components. Similar 


profiles have been observed in the CS survey of NGC1333 by 


Langer et al 


(19961). 


Table El reports the velocity-integrated intensity in each position (the area under the 
Gaussian profile and integrated over the velocity interval at the “zero-intensity” level) 
which we will use to derive the gas temperature, density and species column densities. 
Note that we fitted separately the low (LVG) and high (HVC) velocity components. While 
this is easy for the West lobe emission where the two components are well separated, the 
fit is more tricky in the East lobe, where the two components are blended. In order 
to compute the area of the two components in the East lobe, we precede as follows. 
The LVC component is supposed to be centered on the systemic velocity of the cloud 
(7.6 km s“^) and to have a fixed width of 1 km s“^ (same width as the LVC component 
of the West lobe). The residual spectrum, after removing the LCV component, is HVC. 
Using this method, the resulting LVC integrated intensity is less than 30% of the HVC 
one. Therefore, the uncertainty on HVC intensity produced by the mixing of the two 
components is low whereas the LVC one is less robust. For this reason, we will consider 
in the following only the HVC component of the East lobe. 

Some but not all of the spectra of the West lobe can be fitted by two Gaussians, 
representing the LVC and HVC components respectively. On the contrary, none of the 
spectra of the East lobe can be analysed by two Gaussians. In those cases where the 
two Gaussian fit the line, it would be more correct to consider the line intensity of 
the LVC and HVC components computed in this way, rather than iust integrating over 
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Table 2. Intensity (Tmb) in K, velocity width {Av) in km s“^ and integrated intensity 
(/ TMB(5t>) in K km s“^ of the SO and SO 2 line profiles shown in Fig. [21 LVC and HVC 
mean the low and the high velocity component of the spectra of the west lobe. The 
linewidth of E(LVC) is fixed to 1 km s“^ based on the linewidths of the West LVC 
components because the LVC and HVC components of the East lobe are mixed. 




IRAS2 

E(LVC) 

E(HVC) 

W1 (LVC) 

W1 (HVC) 

W2 (LVC) 

W2 (HVC) 


Aa, A5 (") 

0, 0 

72,-12 

72, -12 

-60, 12 

-60, 12 

-96, 24 

-96, 24 

SO 32 ^ 2i 

Tmb 

2.4T0.09 

2.39±0.07 

2.02T0.07 

2.7±0.08 

1.4T0.08 

2.0±0.08 

0.9T0.08 


Av 

1.1T0.03 

1 

5.0±0.5 

1.1T0.09 

3.2T0.2 

1.1T0.09 

7.2T0.2 


fT mb^v 

3.5±0.6 

3.6±0.6 

14.3±2.3 

3.1±0.5 

6.1T1.1 

1.9±0.3 

7.1T1.2 

SO 23 ^ I2 

Tmb 

<0.3 

0.3±0.07 

0.24±0.07 

0.31T0.05 

0.20±0.05 

<0.15 

0.15±0.05 


Av 

- 

1 

7.7T1.3 

0.8T1.9 

2.8±0.6 

- 

8.3±1.6 


fT mb^v 

<0.6 

0.4T0.08 

2.0±0.7 

0.3±0.1 

0.6T0.2 

<0.3 

1.2T0.4 

so 43 —^ 82 

Tmb 

l.SiO.l 

2.2T0.1 

1.8T0.1 

2.6±0.2 

1.5±0.2 

1.7±0.2 

0.8±0.2 


Av 

1.04T0.06 

1 

6.5±0.2 

0.89±0.06 

2.7±0.3 

0.8±0.1 

7.6T0.3 


fT mbSv 

1.9T0.4 

2.3±0.4 

11.3±1.9 

2.7±0.5 

5.6±1.0 

2.0T0.4 

6.4T1.2 

so 65 ^ 54 

Tmb 

0.69T0.08 

1.7±0.1 

1.2±0.1 

0.56±0.05 

0.77±0.05 

0.24T0.06 

0.38±0.06 


Av 

1.5±0.1 

1 

7.1±0.2 

0.8±0.1 

2.7±0.2 

1 

8.9T0.4 


fT mb^v 

1.2±0.3 

2.4±0.3 

8.3±1.4 

0.44±0.09 

2.6±0.5 

0.26±0.09 

3.6±0.7 

SO2 3i,3 ^ 2o,2 

Tmb 

0.09±0.02 

0.4T0.02 

0.44T0.02 

0.13T0.02 

0.17T0.02 

0.08±0.02 

0.09T0.02 


Av 

2.4T0.7 

1 

4.6T0.4 

l.OiO.l 

4.3±0.3 

l.OTO.l 

9.7±1.1 


fT mb^v 

0.23±0.09 

0.4±0.08 

2.8±0.5 

0.14±0.03 

0.8±0.2 

0.08±0.03 

1.0±0.2 

SO2 10l,9 ^ lOo.lO 

Tmb 

<0.06 

0.08±0.02 

0.18±0.02 

<0.05 

0.06±0.02 

<0.6 

<0.6 


Av 

- 

1 

5.4±0.3 

- 

5.2±0.9 

- 

- 


fT mbSv 

<0.2 

0.1T0.05 

2.8±0.5 

<0.07 

0.25T0.09 

<0.06 

<0.2 

SO2 5i,5 ^ 4o,4 

Tmb 

- 

- 

0.34T0.02 

<0.1 

0.23±0.04 

<0.07 

0.11±0.02 


Av 

- 

- 

5.1±0.4 

- 

2.6±0.4 

- 

12.3±1.0 


fT mbSv 

- 

- 

2.2±0.4 

<0.2 

0.7±0.2 

<0.08 

1.3T0.3 

SO2 62.4 ^ 4i,3 

Tmb 

- 

- 

0.23±0.02 

<0.09 

<0.09 

<0.07 

<0.07 


Av 

- 

- 

6.2T0.3 

- 

- 

- 

- 


f TmbSv 

- 

- 

1.2T0.2 

<0.1 

<0.2 

<0.07 

<0.3 


analysis is not possible everywhere, we decided to adopt the same method for all the 
points. While this choice does not practically affect the integrated intensity of the HVC 
component (the errors remain within the observational errors), it affects the LVC line 
intensity determination. To give the order of magnitude of the possible error associated 
with the adopted procedure, the two Gaussians analysis would give at most a factor 
of 2 lower integrated intensity for the SO lines, 20% less for the CS lines and no LVC 
component would be detected in the SO 2 lines. The consequencies for the computed LVC 
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Table 3. Column densities of SO, SO 2 and CS, as derived from an LVG analysis, at 
several points along the outflow and on the source IRAS2 position The column densities 
(averaged on the 24" beam) have been derived with a non-LTE LVG analysis (see text). 
As we only have one SO 2 transition detected on the source position, we used the same 
temperature and density as for CS to derive the SO 2 column density on IRAS2. LVC and 
HVC refer to the low and high velocity component respectively (see Sect.0]for details). 
The uncertainties in the derived column densities are given by the observational errors, 
and are about 30%. Consequently, the uncertainty in the abundance ratios is about 60%. 
If we compute the LVC SO column densities with the LVG model using a temperature 
of 20 K and a density of 5 x 10^ cm“^ instead of the method explained in section 5.2, 
we find column densities two times higher than the ones reported in this table. 


Aq, A(5 (”) 

Nso 

(lO^"*) 

Nso2 

(10^®) 

Ncs 

(10^®) 

SO 2 /SO 

CS/SO 

0 0 (IRAS2) 

0.8 

2.1 

4.0 

0.2 

0.5 

48 -12 (HVC) 

1.5 

2.5 

0.4 

0.2 

0.03 

60 -12 (HVC) 

3.0 

7.0 

1.2 

0.2 

0.04 

72 -12 (HVC) 

3.5 

7.0 

2.6 

0.2 

0.07 

84 -12 (HVC) 

2.0 

2.0 

1.8 

0.1 

0.09 

-36 12 (LVC) 

0.9 

< 0.5 

4.5 

< 0.05 

0.5 

-36 12 (HVC) 

1.0 

2.0 

3.0 

0.2 

0.3 

-48 12 (LVC) 

0.6 

1.2 

3.8 

0.2 

0.6 

-48 12 (HVC) 

1.0 

2.0 

2.0 

0.2 

0.2 

-60 12 (LVC) 

0.9 

1.2 

3.7 

0.1 

0.4 

-60 12 (HVC) 

1.5 

2.5 

2.5 

0.2 

0.2 

-84 24 (LVC) 

0.6 

0.8 

2.5 

0.1 

0.4 

-84 24 (HVC) 

1.5 

2.5 

4.3 

0.2 

0.3 

-96 24 (LVC) 

0.7 

0.7 

2.7 

0.1 

0.3 

-96 24 (HVC) 

2.0 

6.5 

2.9 

0.3 

0.1 


section and reported in Tables lA. II and nil are the following: Nso would be 30% to 40% 
lower, Ncs would be 20% lower and the SO 2 /SO abundance ratios would have an upper 
limit around 0.08 in the West LVC whereas the CS/SO abundance ratios would be almost 
the same as in Table 01 For this reason, our analysis will especially focus on the HVC 
emission, where the determination of the integrated emission is robust. On the contary, 
we will treat the LVC emission with more caution in the rest of the paper. 


5. Column densities and abundance ratios 


We estimated the column densities of the SO, SO 2 and CS molecules and put some 
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Fig. 4. SO 2 /SO and CS/SO abundance ratios (TableOJ as a function of the distance to 
the protostar. Upper panel represents the LVC ratios of the West lobe (SO 2 /SO ratios 
are upper limits). Lower panel shows the ratios in the HVC of East and West lobes. Note 
that the uncertainty in the values of the ratios is about 60%. 


along the outflow. Both the rotational diagram method and a non-LTE LVC code were 
used. In this section, we present only a summary of the results. The details of this analysis 
are given in Appendix IXI 

For the LVC component, our modeling suggests densities between 2 x 10^ and 2 x 
10® cm“® and temperatures about 25 K. Thus for the modeling of the chemistry (next 
section), we will assume a density of 5 x 10^ cm“® and a temperature of 20 K for this 
component. For the HVC component, we found a difference in the physical conditions 
between the two lobes of the outflow: the temperatures and densities derived from SO 2 
in the East lobe suggest a gradient along the outflow. We also found a slight difference 
between SO and SO 2 . Indeed, SO seems to probe a gas less dense that the gas probed 
by the SO 2 transitions, suggesting a spatial differentiation in the formation of these two 
molecules. However, the data are not sufficient to better quantify this, so that in the 
following we will assume in our chemical modeling that both molecules originate in the 
same gas, where the density is around 10® cm“® and the temperature is about 50 K, the 
median values between SO and SO 2 derivations. 

The column densities derived using both methods differ by less than 30%. The column 
densities and the SO 2 /SO and CS/SO abundance ratios obtained with the LVG model 
are summarized in Tablel^for the several positions in the outflow. The abundance ratios 
are shown as a function of the distance to the protostar in Fig. 01 Remarkably, the derived 
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Table 4. Initial conditions for the chemical modelling. “ Elemental abundances used to 
compute the gas phase composition in molecular clouds (step 1). ^ Species contained in 
grain mantles and sputtered in the gas phase in shocks (step 2). We took the abundances 
observed o n grains in high m ass protostars regions. The cosmic abundance of sulphur 


is 3 X 10® I Sofia et al. 


1994). We used a 30 times lower abundance to take into account 


the hig h depletion of sulphur in molecular clouds IjTieftruuk et al 


imi 


Wa,ke1a,m et a,l 


Ref: (1) 

Wilson & RoodI 1 

1994) 

; (2) 

Cardelli et al. 

(1996 

); (3) 

Mever et 

1. ( 

1998 

); (4) 

Schutte et a 

1, Il996l: (5) 

^eane et al 

J2001r): Ibllchiaret al. 1 

1996 

); (7) I 

looffert et al. 

Jl998ll: (8) I 

^ahimbo et al 

(199 

7J; (9) upper limit from 

van Dishoeck & Blake 

n99sli. 


Elemental abundances^ 


Species 

Abundance 

Ref. 

He 

0.28 

(1) 

0 

6.5 X 10"’* 

(2) 

c+ 

2.8 X 10"'^ 

( 3 ) 

g+c 

1.0 X 10"® 


Grain mantle composition*’ 

Species 

Abundance 

Ref. 

H2O 

10"'^ 

( 4 ) 

H2CO 

4 X 10“® 

( 5 ) 

CH3OH 

4 X 10“® 

(6) 

CH4 

10 "® 

( 7 ) 

OCS 

10 -" 

(8) 

H2S 

10 -" 

( 9 ) 

Sd 

3 X 10"^ - 3 X 10“® 



not only in the kinematics and gas temperature and density, but also in the chemical 
composition. 

In the next section we will try to understand the differences in the chemical com¬ 
position of the LVC and HVC components, as well as the gradients obs erved along the 


outflow, comparing the observed ratios with the theoretical ratios of the 


Wa,ke1a,m et a,l 


1 2004a|l chemical model, predicted for the physical conditions described above. 


6. Chemical model 


In this section, we compare the abundance ratios found in the LVC and HVC components 
with t he predictions the chem ical model we developed recently to study the sulphur chem¬ 


istry llWakelam et al 


2004all . Briefly, the model computes the evolution of the chemical 


composition of a gas at a fixed temperature and density, given an initial composition of 
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In this sense, it is a pseudo time-dependent model, for it does not take into account the 
physical evolution of the gas, but only the chemical evolution. For the specific case of 
the molecular shocks we model here, the shock is assumed to warm and compress the 
gas, giving rise to a sudden increase of the gas temperature and density. Also, the grain 
mantles are sputtered by the ions accelerated in the shock (e. g. Pineau de Foret s et al. 
1993) and/or destroyed by the collisions between grains (e.g. 


Jones et al. 


1994 . As a 


result, some or all the grain mantle components are injected into the gas phase at the 
shock passage. The model thus computes the gas chemical composition as a function of 
time, following those changes: gas density and temperature, as well as the gas chemical 
composition due to the grain mantle (partial) evaporation. 


Note that the model does not aim to reproduce all the molecular abundances, but it 
specifically focuses on the sulphur chemistry. In this respect, it has the most up-to date 
and complete set of reactions involving S-bearing molecules. 


In practice, in the present study the passage of the shock is simulated by a two-step 
process: 

1) Before the shock, the gas composition is computed for a cloud with the density and 
temperature derived by the observations of the LVC, namely 5 x 10^ cm“^ and 20 K. 
The elemental abundances used for this step are summarized in Table 01 The age of the 
cloud is estimated by comparing the model predictions with the observed CS/SO and 
SO 2 /SO abundance ratios. 

2) The cloud chemical composition computed in step 1 is used as the initial condition 
for step 2 computations. To the cloud composition, we add the grain mantle components 
released into the gas phase (see Table 0)1, and then the model follows the chemical 
evolution of the gas at 0.1 —10 x 10® cm“® and 25-120 K (namely a range of temperatures 
a nd densities app r opriate for the conditions found in the HVC). Following the results 

Wakelam et all lj20Q44 . we assume that sulphur mainly evaporates in the H 2 S, OCS 


of 


and atomic form and we vary the amount of atomic sulphur released in the gas phase. 
Note that we also considered the possibilit y that sulphur evapora tes in the molecular 


form S 2 but as in the case of the hot cores llWa.kela.m et af 


20043 ) the model with the 


evaporation of S 2 does not reproduce the observations in shocks. We will therefore not 
discuss this case further. 


In summary, we studied the influence of four parameters: the temperature, the density, 
the time since grain mantle sputtering and the amount of atomic sulphur released from 
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Fig. 5. Theoretical SO 2 /SO and CS/SO abundance ratios as a function of time. The 
temperature is 20 K, and the density 5 x 10^ cm“^. The two grey bands represent the 
observed ratios in the LVC component of the East lobe (maximum and minimum values 
from Tablel^with the uncertainty of 50%). Note that for the SO 2 /SO ratio, we only have 
a upper limit in the -36", 12" LVC position. The darker band corresponds to the overlap 
of the observed values. 

6.1. Results: the cloud gas (LVC component) 

Fig. 0 shows the SO 2 /SO and CS/SO abundance ratios as predicted by the chemical 
model for the LVC temperature and density (20 K and 5 x 10^ cm“^) as a function 
of time. The observed ratios are reported on the same figure. The comparison between 
the model predictions and the observed ratios indicate that the cloud age is around 
2 X 10® yr or greater than about 10® yr. In that case, the observed SO 2 /SO ratios are 
in agreement with the model whereas the predicted CS/SO ratios are four times lower 
than the observed ones. An increase of the density or a decrease of the temperature 
leads to a better agreement between observed and predicted ratios. A decrease of the 
density or an increase of the temperature would have the opposite effect. For an age of 
2 X 10® yr, the model predicts absolute abundances 10 times higher than the observed 
ones whereas the observed abundances are reproduced by the model within a factor of 3 
for times longer than 10® yr. This analysis assumes that SO, SO 2 and CS all originate in 
the same gas, which in the case of the molecular cloud is likely a good approximation. 
As shown in Fig. [3 the CS/SO and SO 2 /SO ratios are each clocks that work differently 
depending on the temperature and density. Indeed, because of the different dependence 
on the different model parameters, the use of both ratios can help constrain the age along 
with the density and temperature of the gas, as discussed above. 

6.2. Results: the shocked gas 

As we already mentioned, in addition to the temperature, the density and the time, 
we studied a fourth parameter: the amount of atomic sulphur injected in the gas phase 
which depends on the shock strengh. This last parameter, the amount of atomic sul- 
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[S]ini 

Fig. 6. Theoretical SO 2 /SO (solid line) and CS/SO (dash dot line) abundance ratios 
as a function of the initial amount of atomic sulphur evaporated in the gas phase. We 
varied the amount of atomic sulphur (abundance compared to H 2 ) evaporated in the 
gas phase between 3 x 10“^ and 3 x 10“®. All the other model parameters are fixed: 
T=50 K, n(H2)=10® cm“^ t=2000 yr. The two grey bands represent the observed ratios 
in the shocked gas of the East and West lobes (see Table El). Note that the two bands 
are superimposed in the middle (Log(Ratio)=0.004-0.1). 


dependence of the predicted abundance ratios on this parameter. We took the average 
density and temperature previously derived by the SO emission in the HVC, namely 
10® cm“® and 50 K respectively. Anticipating the results on the time dependence of the 
predicted ratios, we took an early time, 2000 yr. For these conditions, we ran several 
models, varying the initial amount of atomic sulphur. Fig. El shows the results of the 
modelling together with the range of observed ratios in the two lobes of the outflow. The 
first robust result is that only a small fraction of the sulphur is evaporated in the gas 
phase: < 1/100 of the elemental sulphur in the West lobe and < 1/40 of the elemental 
sulphur in the East lobe. This suggests that the shocks in IRAS2 are rather mild. The 
second result is that H 2 S and OCS evaporate from the grain more easily than atomic 
sulphur. Indeed, decreasing the OCS injected in the gas by the same amount as S would 
give predictions not in agreement with the observations, because the abundance of CS 
in the warm shocked gas highly depends on the initial abundance of OCS. The decrease 
of H 2 S does not influence the results signihcantly. We will discuss this point further in 
Sect. 7. 

Taking the above composition for the injected material, we explored the dependence 
of the SO 2 /SO and CS/SO abundance ratios as function of the other three parameters 
of the model: the density, temperature and age of the shock. The results are shown in 
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Fig. 7. Theoretical SO 2 /SO and CS/SO abundance ratios as a function of three param¬ 
eters: time (upper panel), the H 2 density (middle panel) and the temperature (lower 
panel). In each panel, the other parameters are hxed: T=50 K and n(H2)=105 cm“^ in 
the upper panel, T=50 K and t=2000 yr in the middle panel, and n(H2)=105 cm“^ and 
t=2000 yr in the lower panel. For all the figures, the initial amount of evaporated atomic 
sulphur is 3 X 10“^. The two grey bands represent the observed ratios in the shocked gas 
of the East and West lobes, (see Table 0). Note that the two bands are superimposed in 
the middle (Log(Ratio)=0.004-0.1). 

The SO 2 /SO ratio depends weakly on the shock age, until 10^ yr, and on the density 
and temperature of the shocked gas. On the contrary, the CS/SO ratio increases with time 
for (2 — 3 ) X 10 ^ yr before decreasing sharply because SO is efficiently produced whereas 
the abundance of CS does not change. The CS/SO also decreases with the density, for 
’i^H 2 ^ 10^ cm“^ and with the temperature. The CS/SO ratio is therefore an efficient 
constraint on the shock age, temperature and density. We will discuss the constraints 
derived along the outflow in Sect. 7. 

Finally, we notice that the HVC chemical modelling depends on the initial gas phase 
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too low to reset the chemistry as is the case for the hot cores. In addition, the CS/SO 
ratio depends more on the initial composition than the SO 2 /SO ratio because CS is also 
directly influenced by the initial abundance of atomic carbon (wich is a function of the 
cloud age). Indeed, the initial abundances are parameters hard to constrain, which we 
decided not to alter. The absolute abundances depend even more on the model parameters 
than the abundance ratios. Indeed, we are not able to reproduce the observed absolute 
abundances which are two orders of magnitude higher than the modelled ones. This 
discrepancy can have various causes. The first is that the modelling is wrong, and only 
further observations will be able to rule out this possibility. Other possibilities include the 
uncertainty on the adopted density and temperature. Here we are using average values 
across each modeled point, but there may be a stratification of density, temperature and 
chemical composition, that may affect the results. Also, the derived absolute abundances 
are rather uncertain because they are based on observations of low J CO transitions 
and this may be the most important cause of the discrepancy between the model and 
the observations. These observations have relatively large beams, and the amount of gas 
probed by the CO lines can be substantially larger than the gas where the SO, SO 2 and 
CS lines are emitted, even up to a factor of 100. Only high spatial resolution observations 
will be able to clartify this possibility. In addition, the CO 2-1 line used to derive the H 2 
column density may be optically thick, underestimating the H 2 column density, and only 
^^CO observations will be able to clarify this point. In the absence of such observations 
we think that our modelling, based on the molecular abundance ratios -which are more 
robust than the absolute abundances- gives a reasonable and consistent interpretation of 
the observations. We will try to explore in detail the consequences of this interpretation 
in the next section. 


7. Discussion 

7.1. Nature of the LVC component 


The density and temperature derived in the LVC (Sect. roil are consistent with the 
hypothesis suggested in §4, that the LVC is associated with the molecular cloud. Indeed, 
observations of the extended ^^CO and ^^CO em ission suggest a density of ~ 10^ cm~ ^ 
and temperature of ~ 20 K for the gas cloud (see 


Warin et al. 


1996; 


Bergin et al 


2003). 


Likely, this gas is rather in the warm layers just behind the PhotoDissociatio n Region ere 


Bergin et a,l 


ated by the illumination from a Go ~ 100 — 400 FUV held (see discussion in 
2i]02|). The SO column density in the LVC component is about 5 — 9 x 10^^ cm“^, which 


ers extend for about 8 x 10^^ cm ^ as derived by 


6-1x10-9, a 

ssnming that the warm lav- 

Warin et al. 

(1996 

). This SO abundance 
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19951) and the modeling of them I Wakelam et al.M2004a^ . so that we are tempted to 
conclude that the LVC component is dominated by the cloud emission. 

Also from a chemical point of view, the SO 2 /SO and CS/SO ratios observed in LVC 
seem to reflect the cloud chemistry, even though the chemical model predicts a CS/SO 
ratio lower than what is observed. 

Finally, the constraints on the cloud age given by the sulphur chemistry , i.e. ~ 2 x 


10 ^ yr or > 10® yr, are remarkably in agreement with previous estimates from 


Lada et al 


(199^, < (1 — 2) X 10® yr, obtained by analysing Near-IR survey data. Therefore, the 
cloud hosting 1RAS2/NGC1333 is chemically young, less than 10® yr old. 


7.2. Age and physical conditions of the outflow 


The comparison between the observed abundance ratios and the chemical modelling 
presented in Sect. 6 gives some constraints on the shock age, density and temperature. 
First, the sulphur bearing species indicate an age of 5 x 10^ — 7 x 10® yr for the outflow 


Bachiller et al 


1199 ^, 


of 1RAS2. This age is in agreement with the range estimated by 
4 X 10^ — 5 X 10® yr, based on the dynamics of the bullets of the 1RAS2 outflow. The 


estimated age of the outflow is also in agreement with the genera l schem e of formation 


and destruction of molecules in shocks proposed by 


Wakelam et al 


1 2004ai , Table 7): the 


simultaneous presence of SO and SO 2 would indicate a relatively young shock (< 10® yr). 
The accuracy on the SO 2 /SO and CS/SO ratios is however not good enough to constrain 
the age of each clump and therefore to test if they are different events. The density of the 
HVC component in the West lobe is constrained by the CS/SO ratios to be lower than 
2 X 10® cm“® and the temperature lower than 70 K. These physical conditions are in 
remarkable agreement with the non-LTE analysis of the SO emission presented in Sect. 5 
(Table ESI). No constraints on the temperature or density can be deduced for the HVC 
component in the East lobe. However, the chemical modelling shows that an increase of 
the density and/or of the temperature leads to smaller CS/SO ratios. Thus, the lower 
CS/SO ratios observed in the East lobe compared with the West part of the outflow 
can be explained by a denser and/or warmer gas. The study of molecular excitation 
conditions (Table ESI seems to confirm that the gas in the East lobe is denser than in 
the West lobe. 

The density and te mperature we derived in HVC are lower than that obtained by 


Jorgensen et al 


2004 70 K and 10® cm ®). For the temperature, we are still within the 


mutual uncertainties. For the density value, we are forced to conclude that 10® cm“® is not 
in agreement with the HVC SO and SO 2 emission data since a high density would require 
gas temperatures too low, similar to the rotational temperatures (in Table IXtIi . Indeed, 
we showed that optically thick lines cannot explain these low rotational temperatures (see 
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that they did not consider the two lobes of the outflow separately, and also they treated 
all the studied molecules without studying the excitation conditions of each molecule 
individually. As we showed, the physical conditions are different in the East and West 
lobes, and different molecules probably trace different layers of shock ed gas (see also 


Bacdnller et al.l200lll . In any case, the SO and CS densities computed by 


Jorgensen et al 


are similar to ours in both LVC and HVC. 

The low temperature (< 70 K) derived in t he HVC compon e nt fro m SO and SO 2 


emission, and from other molecular emission by 


Jorgensen et al 


1 200411 ■ as well as the 


small fraction of sulphur released in the shocked gas found from the chemical modelling 
suggest mild shocks. The SO, SO 2 and CS molecules seem rather to trace the entrained 
material, mildly shocked, rather than the shock from the outflowing wind. This would 
explain both the lack of gas compression and the low gas temperature. 


7.3. Consequences for the depleted form of sulphur 


In this work on sulphur chemistry in shocks, we assumed that sulphur was sputtered 
from the grains in the OCS, H 2 S and S forms. This assumption does not imply that 
sulphur is depleted in the atomic form but that the third form of depleted sulphur is 
quickly convert e d into S once evaporated in the gas phase. We denote this species S?.In 


Wakelam et al 


l 2004a|l . we proposed that S? is a polymeric form of sulphur such as Sg. 


One conclusion of this work is that it is more difficult to sputter S? than OCS and 
H 2 S. The efficiency of each molecule sputtering depends on the binding energy of that 
molecule. Thus S? should have a higher binding energy than OCS and H 2 S. Unfortunatly, 
this parameter is poorly known and depends highly on the grain surface composition. 


8. Conclusions 

We studied the SO, SO 2 and CS emission in the protostar IRAS2 region and its East- 
West outflow. Using both LTE and non-LTE methods, we computed the species column 
densities of the gas associated with the outflow and the ambient medium respectively. 
The computed CS/SO and SO 2 /SO abundance ratios were compared with the predictions 
of a chemical model. 

The results are: 

- The sulphur bearing emission is similar in the IRAS2 position and in the ambient 
medium whereas it is enhanced in the shocked regions at the interface between the 
cloud and the wind. 

- The low velocity component (LVC) of the sulphur bearing spectra reflects the physical 
conditions and the chemistry of the cloud younger than 10® yr. 

- The chemical modelling suggests that the IRAS2 outflow is young (< 2 x 10® yr) but 
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time each clump of the outflow. The density of the gas traced by SO and SO 2 is lower 
than 2 x 10® cm“® and the temperature lower than 70 K in the West lobe. The gas 
density in the East lobe is higher. 

- The theoretical modeling suggests than SO and SO 2 trace different layers of gas 
and the entrained material rather than the shock itself. Thus, these molecules are 
not appropriate to study the physical conditions of the gas and the evolution of the 
outflow. CS on the contrary seems to be more sensitive to these parameters and the 
use of CS/SO ratios can help in understanding the molecular shocks associated with 
the outflows. 

- Sulphur-bearing species depleted on the grain mantles may be more difficult to sputter 
than OCS and H 2 S. Thus, this species should have a higher binding energy than OCS 
and H 2 S. 
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Table A.l. Rotational temperatures and column densities of SO and SO 2 computed 
from the rotational diagrams of Fig. ro Note that the detected lines of SO 2 transitions 
towards IRAS2, Wl(LVC) and W2(LVC) are not numerous enough to construct the SO 2 
rotational diagram in these regions. Thus the SO 2 column densities towards IRAS2 and 
the LVC components (W1 and W2 LVC) are derived assuming LTE and an excitation 
temperature of 10 K. 


Aa, A5 (") 

Trot{K) 

SO 

Aso(10i%m-2) 

TrotiK) 

SO 2 

AsoJlO^'cm-^) 

CO 

Aco(10i%m-2) 

0, 0 (IRAS2) 

11 ±3 

0.9 ±0.01 

10 

0.5 ±0.2 

- 

72, -12 (E HVC) 

13 ±3 

4.0 ±0.4 

28± 11 

10 ± 1 

0.75 

-60, 12 (W1 LVC) 

7.3 ±0.2 

0.8 ± 0.04 

10 

0.3 ±0.07 

- 

-60, 12 (W1 HVC) 

11.7T2.8 

1.6±0.2 

28.1±12.9 

2.8±0.2 

0.45 

-96, 24 (W2 LVC) 

7.2 ±0.6 

0.6 ±0.07 

10 

0.2 ±0.07 

- 

-96, 24 (W2 HVC) 

12.5±1.6 

2.6±0.1 

19.6±7.1 

3.2±0.3 

- 


MNRAS, 281, L53 

Warin, S., Castets, A., hanger, W. D., Wilson, R. W., & Pagani, L. 1996, A&A, 306, 935 
Wilson, T. L. & Rood, R. 1994, ARA&A, 32, 191 


Appendix A: Determination of the column and physical densities 

In this section we report the details on the estimates of the column densities of the SO, 
SO 2 and CS molecules at selected positions along the ouflow. We used both the rotational 
diagram method and a non-LTE LVG code. The two methods and the relevant results are 
discussed separately in the following. We anticipate that they give values for the column 
densities that differ by less than 30%. 


A.l. Rotational diagram method 


As widely discussed in the literature (e.g. 


Goldsmith & hanger 


199S), the rotational di¬ 


agram method rests on the assumption that the molecular levels are in local thermal 
equilibrium (hTE) and that the lines are optically thin. Using the four observed transi¬ 
tions of SO and SO 2 , we constructed rotational diagrams at the three positions E, W1 
and W2, and at the IRAS2 position (see Fig. Ell). Note that we did not perform this 
analysis for the CS data because we have only two transitions available. The obtained 
SO and SO 2 rotational temperatures and column densities are reported in TablelA.il 


As previously found for CH 3 OH llSandell et al. 


19M 


Bachiller et al 


1998), the ro¬ 


tational temperatu res are relativ ely low when compared to the expected temperatures 


in the shocked gas. 


Sandell et al 


concluded that IRAS2 was the first example of a cold 
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Eupp/k (K) 


Fig. A.l. Rotational diagrams of the SO and SO 2 molecules in the direction of IRAS2, 
E, W1 and W2 regions. In W1 and W2 we have constructed separate diagrams for the 
LVC and HVC components: Wl(LVC), Wl(HVC), W2(LVC) and W2(HVC) whereas in 
E, we give only the HVC component (see text). Stars and triangles correspond to the SO 
and SO 2 observed data respectively. 


Indeed, the gas temperature in the IRAS2 outflow is certainly higher (i.e. > 60 K) 


detected towar ds the two lobes (see 


terpretation by 


Bachiller et al 


TablelA.il since 

bright NH?if3.3l emission has been 

Bachiller et al. 

1998 

). We therefore endorse the in- 


and we demonstrate that this is also the case for the SO 


and SO 2 molecules in Appendix m the low rotational temperatures are simply due to 
non-LTE effects. 


A.2. Non-LTE LVG method 


In order to take into account the expected non-LTE effects, we analysed the data by 
means of a non-LTE LVG (Large Velocity G radient) code. The general details of the 


code are described in 


Ger.carelli et al 


1 2flfl2|l . We included the molecular data and the 
collisional coefficients for the SO, SO 2 an d CS molecules, by taking the former fr om the 


JPL Gatalogue (http://svec.jpLnasa.qov/; 


Pickett et al 


1991 ) and the latter from Green 
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Log(Density [cm ) 


Fig. A.2. contour at 2a (95.4% in confidence) of the LVC and HVC SO line emission 
in the 5 positions of the West lobe (upper panels) and the HVC SO and SO 2 line emission 
in the 4 positions in the East lobe (lower panels) compared to the non-LTE LVG code. 


LVG code computes the rotational line fluxes solving simultaneously the statistical equi¬ 
librium of the levels and the radiative transfer of the emitted photons. The approximation 
of the LVG consists of assuming that the photons emitted at any point of the medium 
either are absorbed in the immediate vicinity or they escape (this is the so-called photon 
escape formalism), because the different regions are moving at high velocities. In prac¬ 
tice, the photons cannot be absorbed if their frequency -shifted by the Doppler effect- 
is outside the thermal width centered on the rest frequency of the transition. The LVG 
approximation assumes homogeneous and isothermal semi-infinite slab. Although this 
approximation is strictly correct in the presence of large velocity gradient, and the ab¬ 
sence of temperature and density gradients, it is very often used to give a first estimate 
of the average gas density and temperature of the observed medium. In this work, we 
will use the non-LTE code in this way. More sophisticated models would be required for 
a more sophisticated analysis. However, given the quality of the obtained data, a greater 
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Table A.2. Temperatures and densities constrained by the SO and SO 2 observed line 
emission using a non-LTE LVG method (see text) for the different positions along the 
outflow. The last two columns summarise the temperatures and densities used to compute 
the CS column densities reported in Table 0 


Aa, A5 (") 

T(K) 

SO 

n (cm“®) 

T(K) 

SO 2 

n (cm“®) 

Adopted for CS 

T (K) n (cm-®) 

0 0 (IRAS2) 

- 

> 2 X 10"^ 

- 

- 

20 

5 X 10* 

48 -12 (HVC) 

- 

> 2.5 X 10"^ 

< 35 

> 6.3 X 10® 

50 

10® 

60 -12 (HVC) 

> 30 

3.1 X 10"* - 4.0 X 10® 

20-80 

2.5 X 10® - 6.3 X 10® 

50 

10® 

72 -12 (HVC) 

- 

> 5 X 10"* 

20-50 

5 X 10® - 7.9 X 10® 

50 

10® 

84 -12 (HVC) 

- 

> 2.5 X 10® 

< 30 

> 1.2 X 10® 

50 

10® 

-36 12 (LVC) 

14 - 75 

2.2 X 10"^ - 2.5 X 10® 

- 

- 

20 

5 X 10* 

-36 12 (HVC) 

- 

> 2.5 X 10‘‘ 

- 

- 

50 

10® 

-48 12 (LVC) 

14 - 130 

2.2 X 10'* - 5.0 X 10® 

- 

- 

20 

5 X 10* 

-48 12 (HVC) 

> 70 

4 X 10* - 1.6 X 10® 

- 

- 

50 

10® 

-60 12 (LVC) 

15 - 24 

1.0 X 10® - 2.2 X 10® 

- 

- 

20 

5 X 10* 


> 25 

2.8 X 10* - 6.3 X 10* 





-60 12 (HVC) 

> 45 

2.5 X 10* - 1.2 X 10® 

- 

- 

50 

10® 

-84 24 (LVC) 

< 130 

> 2.2 X 10* 

- 

- 

20 

5 X 10* 

-84 24 (HVC) 

> 60 

2.8 X 10* - 1.0 X 10® 

- 

- 

50 

10® 

-96 24 (LVC) 

< 55 

> 2.5 X 10* 

- 

- 

20 

5 X 10* 

-96 24 (HVC) 

- 

> 2.5 X 10* 

- 

- 

50 

10® 


The code has three free parameters to be constrained by comparing the model pre¬ 
dictions with observed data: the gas density and temperature and the column density of 
the species. Implicitly, we are assuming that the emission fills the beam (24") and the 
observed linewidth, namely 4 and 5 km s“^ for HVC of SO and SO 2 and 1 km s”’^ for 
LVC of both molecules. We ran several cases varying the SO and SO 2 column densities 
from 10^^ to 10^^ cm“^, the density between 10"^ and 10^ cm“^ and the temperature 
between 10 and 200 K. We then searched for the minimum in this 3-D space, where 
the is defined as usual : 

2 1 {Observations — Model)'^ 

^ ^ N-2 ^ ^ 

We computed the SO and SO 2 column densities, as well as the temperature and 
density of the gas in five and four positions of the west and east lobes of the outflow 
respectively (Tables 13 and rot . These positions include the E, W1 and W2 regions 
discussed in previous sections. 

Taking the species column densities with the lowest absolute (Table 0) , we con¬ 

strained the gas density and temperature by plotting the contour at 2a (95.4% in 
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temperature and density can reproduce the observed fluxes, so that they are only ap- 
proximatively constrained. Tablereports the kind of constraints we obtained at each 
studied position. 


Having available only two transitions for CS we could not carry out a similar method 
to constrain the gas density and temperature and the CS column density simultaneously. 
We assumed instead the average density and temperature as suggested by the SO ob¬ 
servations and reported in Table roi (see the detailed discussion in section 5.3), and 
computed the CS column density consequently using the LVG model. The results are 
reported in Table |21 Note that CS having the same density and temperature as SO im¬ 
plies that the two species are also chemically related. This is a major assumption, which 
is partially supported by the theoretical predictions that, in dense an d warm gas, CS 
is mainly formed by C -I- SO (thus SO and CS are chemically linked; 


Wakelam et ai 


Cha.rnlev 




2004all . However, even our observations cast some doubt on the robust¬ 


ness of this assumption, for the SO and CS emissions in the maps are slightly displaced. 
Not having other possibilities, we will adopt this assumption, but in the following we will 
keep in mind the limits of the assumption used. 


A.3. Modeling results 

Gas temperature and density 

The non-LTE analysis gives an indication of the gas density and temperature along 
the outflow and in the LVC and HVC components (Fig. lA.2l and Table roi . within the 
limits discussed in the previous section and summarized again at the end of this section. 
In the following, we discuss the two velocity components separately. 


i) LVC component 

Although it is difficult to constraint exactly (with the presented observations) the gas den¬ 
sity and temperature in this component, the plot in Fig. roi suggests densities between 
2 X 10"^ and 2 x 10® cm“® and temperatures somewhat higher than 15 K. Complementary 
CO J=2-l line observations of the IRAS2 outflow indicate brightness temperatures of 
20-22 K, which implies a gas kinetic temper ature of about 25 K (Wakelam 2005 in prep). 

(199^ who derived a similar temperature (about 


This is consistent with 


Warin et al 


20 K) and a similar density (about 10® cm“®) for the ambient gas by using ^®CO and 
C^®0 observations with a larger beamsize (2.5 arcmin). Thus, for the modelling of the 
chemistry (next section), we will assume a density of 5 x 10"^ cm“® and a temperature of 
20 K for this component. 


a) HVC component 

Here, the temperatures and densities are higher than in the LVC component and different 
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lobes (East and West), and as derived from SO and SO 2 respectively, present a slightly 
different behaviour. The temperature and density derived from SO 2 in the East lobe 
suggest a gradient along the outflow. The two positions of maximum SO 2 emission at 
(60", —12") and (72", —12") display the highest kinetic temperatures, whereas the other 
points, at the edge of the emission peak, display higher densities and lower temperatures. 
The gas probed by the SO 2 emission seems therefore warmed at the peak position, and 
compressed around it. A similar behavior may also be present in the gas probed by the 
SO transitions but it is less evident. Remarkably, Fig. roi shows that the T-n curves 
of the SO transitions do not overlap the curves of the SO 2 transitions. The SO data 
seem to probe a gas less dense than the gas probed by the SO 2 transitions, suggesting 
a spatial differentiation in the formation of these two molecules. Indeed, the map of the 
emission extent of the HVC component suggests that the SO emission is more extended 
than the SO 2 . However, the data are not sufficient to better quantify this, so that in the 
following we will assume in our chemical modeling that both molecules originate in the 
same gas, where the density is around 1 x 10^ cm“^ and the temperature is about 50 K, 
the median values between the SO and SO 2 derivations. The same applies for CS, which 
is also assumed to originate in the same gas the SO and SO 2 . As previously mentioned 
(§5.2), this is likely a rough approximation, but the available data do not allow a better 
refinement of the treatment. 


Hi) Column densities and abundance ratios 

As anticipated, the column densities computed with the non-LTE LVG method (Table|31) 
are very similar to the ones derived by the rotational diagrams. They differ by less than 
30% with the exception of the SO 2 column density in W2(HVC) where the difference 
reaches 50%. The uncertainty on this latter value is higher because two transitions of 
SO 2 were not detected in W2. Note that the SO and SO 2 transitions are predicted to be 
optically thin at the derived column densities. 

Using the derived column densities, we estimated the abundance ratios reported in 
TableOland shown as function of the distance from the protostar in Fig.^] The absolute 
abundances were derived from the CO column densit y, assuming a standard CO to H 2 
ab undance \CO \ = 10 ~^, following 


bvlWarin et al 


Warin et al 


1 1996|1 . We adopted the value estimated 


1 19Q6 i for the LVC : 8x 10^^ cm For the HVC, we took the CO column 


densities obtained from complementary observations of the CO J=2-l line for our study 
of the SiO emission in the IRAS2 protostellar outflow (Wakelam 2005, in prep). The CO 
column density was estimated assuming the levels to be populated according to LTE and 
the transition to be optically thin. This assumption is motivated by the fact that the 
CO profiles are similar to the CS 2-1 ones with a line width of 10-20 km s“^. However, 
only ^^CO 2-1 observations could confirm this hypothesis and the real H 2 column density 
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for the CO gas. The values in the HVC are close to 10^^cm“^; they are summarized in 
Table 4 and are in agreement with the values obtained by Bergin et al. (2003), averaged 
over a region of 4 arcmin size. 

We derive the following average absolute abundances for SO, SO 2 and CS respec¬ 
tively : 9 X 10“^°, 10“^° and 4 x 10“^° in the LVC. We obtain much higher values in the 
HVC : 2 X 10-’’, 4 x 10“® and 4 x 10“® in the West lobe, and 3 x 10"^, 7 x 10“® and 
2 X 10“® in the East lobe. Unfortunatly, we cannot perform an accurate analysis at each 
point because of the lack of observations. Therefore, we will use the abundance ratios for 
the subsequent modelling. 

Remarkably, the derived abundance ratios are different in the LVC and HVC 
components, which, therefore, differ not only in the kinematics and gas temperature 
and density, but also in the chemical composition. The CS/SO ratio is 0.3-0.6 in the 
LVC component, with the ratio steadily decreasing going outward (across the outflow). 
On the contrary, the CS/SO ratio in the HVC component is much more variable: it is 
between 0.03 and 0.1 in the East lobe, and 0.1-0.3 in the West lobe, i.e. up to 6 times 
larger in the West lobe. Contrary to the LVC component, it increases going outward 
in both lobes, with the exception of the farthest point in the West lobe. However, the 
largest difference between the LVC and HVC components is seen in the SO 2 /SO ratio; 
less than 0.005 in the LVC component and about 0.2 in the HVC component, almost 
constant along the outflow. 

iv) Uncertainties 

One has to keep in mind the various approximations/assumptions used when trying 
to interpret these results. First, in our analysis we assumed that the three molecules 
originate in the same gas, which may not be the case. Second, in practice we assumed that 
the emission fills the beam, or that the filling factor is the same for the three molecules. 
If the lines are optically thin this would give values (gas density and temperature, as 
well as the species column density) that are correct, but if the lines are optically thick 
this may not be the case. In practice, we estimate that the computed abundance ratios 
are likely uncertain by about 60%. Higher spatial resolution observations are required for 
better estimates. 


Appendix B: Excitation conditions 


In this appendix, we present a stu dy of the SO and SO-;> e xcitation in the shocked gas using 


a LVG model and following the 


Bachiller et al 


( 1998ri idea. In practice, we would like 


to understand what physical conditions can explain an underestimation of the rotational 
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J4b67ji4b67 
Log(density[cm"-^]) Log (density [err--^J) 


Fig.B.l. SO (left panel) and SO 2 (right panel) rotational temperature as a function of 
density computed for several kinetic temperatures between 60 and 100 K. We assumed 
the SO and SO 2 column densities found in Wl(HVC) (Nso = 1-5 x and 

Ns 02 = 3 X 10^^CTO“^) using rotational diagrams. See text for details. 



Nso [cm"’) L09(derisity [ein"^]) 



Nso: (cm-2) l09(clerslty (em-J)] 

Fig.B.2. Left pannels: SO and SO 2 rotational temperatures as a function of column 
densities computed for several H 2 densities between 10^ and 10^ cm“^and a kinetic tem¬ 
perature of 100 K. Right pannels: Error on the SO and SO 2 column densities computed 
with rotational diagrams as a function of H 2 densities and for several column densities. 

the robustness of the column densities derived with the rotational diagrams if the basic 
hypothesis of the method is not verified. 

Using the SO and SO 2 fluxes computed by the LVG model, we reconstruct theo¬ 
retical rotational diagrams for different column densities, kinetic temperatures and H 2 
densities. From these rotational diagrams, we derived theoretical rotational temperatures 
and column densities. To compare the theoretical diagrams with the observed ones, we 
constructed them only using the LVG fluxes of the four observed transitions of SO and 
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librium (LTE) and of optically thin lines. We studied the effects of the first assumption 
by varying the H 2 density whereas the effects of the optical depth depend on the column 
density. 

Effects of the H 2 density: 

Assume first that the lines are optically thin (Nso and Nso 2 < 10^^ cm“^). For a 
fixed column density, we computed the theoretical Trot varying the kinetic temperature 
between 60 and 100 K and the H 2 density between 10^ and 10^ cm“^ and show them on 
Fig. EH If the kinetic temperature of the gas is higher than 60 K, the computed Trot 
can be as low as the observed ones (7-13 K for SO and 20-30 K for SO 2 , see Table roi 
if the H 2 density is lower than 10® cm“®. The critical densities of the higher transitions 
for the two molecules are around 10® and lO’’ cm“® respectively. So if the gas density is 
lower than 10® cm“®, the LTF is not reached and the higher levels of the transitions are 
not populated, resulting in Trot < Tk. 

Effects of the optical depth: 

To study the effect of the optical depth, we varied the column density between 10^® 
and 10^^ cm“^ (see Fig. ra . We assume LTE (i.e. uhs = lO’^ cm ®). The computed 
rotational temperatures decrease very quickly for SO and SO 2 column densities higher 
than 10^® cm“^. If the kinetic temperature of the gas is higher than 100 K and the LTE 
reached, a SO and SO 2 column density as high as 10®^ cm“^ is needed to obtain Trot as 
iming an H 2 column density in the shocked gas of 5.8 x 10^® cm“^ 
, this gives an SO and SO 2 abundance of 10“"^, which is not realistic. 

Effects on the computed column densities: 

In case of non-LTE and/or optically thick lines, the column densities computed with 
rotational diagrams can be wrong. To quantify the error, we computed the theoretical 
column densities for different real column densities, kinetic temperatures and H 2 den¬ 
sities. In Fig. EH we reported the ratios between the real column densities (used to 
run the LVG model) and the computed ones (derived from the theoretical diagrams) for 
all parameters. Whatever the physical conditions, SO is underestimated. The error is a 
factor of 1.5 if Nsoreai < 10®® cm“^ whatever nua, a factor of 10 and 6 at low and high 
densities respectively if Nsoreai = 10®'’ cm“^. The error is much higher on Nsoa^ for 
Nso 2 reai < 10®® cm“^, Nso 2 computed IS Underestimated by a factor of 6 at low densities 
and overestimated by a factor of 3 at LTE. 

Thus, the most probable explanation for the low observed Trot is a non-LTE effect 
because of H 2 densities lower than 10® cm“®. Considering the observed column densities 
arounf 10®^ cm“® (Table IHJ, the error on the computed column densities should not be 
high: less than a factor of 2 for SO and 6 for SO 2 depending on the H 2 density. 


low as observed. Assi 


( Rachiher et abllOOS 



